First-order Synchronization Transition in Locally Coupled Maps 
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We study a family of diffusively coupled chaotic maps on periodic d-dimensional square lattices. 
Even and odd sub-lattices are updated alternately, introducing an effective delay. As the coupling 
strength is increased, the system undergoes a first order phase transition from a multi-stable to a 
synchronized phase. Further increase in coupling strength shows de-synchronization where the phase 
space splits into two ergodic regions. We argue that the de-synchronization transition is discontin- 
uous for piece-wise linear maps, and is continuous for non-linear maps which are differentiable. 

PACS numbers: 05.45.-a, 05.45.Xt 
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Synchronization is observed in a wide class of complex 
systems. Typically, it appears when the range of the cor- 
relations inside the system is of the same order as the 
system size. Few examples of mutual synchronization in 
complex dynamical systems are flashing fireflies elec- 
tronic circuits and chemical reactions In recent 
years, synchronization of spatially extended systems has 
drawn considerable interest. In particular coupled map 
lattices (CMLs) initially introduced as simple mod- 
els of spatio-temporal chaos, has received a great deal 
of attention as a model of synchronization. It has been 
realized that two different replicas of the same CML are 
coupled directly |^ or through a common external ran- 
dom noise |^ , can become synchronized for large enough 
coupling strengths. Recently, Ahlers and Pikovsky 
pointed out that synchronization transitions in one di- 
mensional CMLs are generically in the universality 
class of the Kardar-Parisi-Zhang (KPZ) model, however 
strongly nonlinear maps could be in the universality class 
of Directed percolation (DP). 

From a point of view of general statistical theory syn- 
chronization is a non-equilibrium phase transition where 
distinct patches of the CMLs oscillate in phase. The 
phenomenon is similar to the roughening of growing in- 
terfaces and can be modeled by multiplicative noise par- 
tial differential equation[^. This picture is further modi- 
fied by Miinoz and Pastor-Satorass|9j to incorporate first 
order phase transitions (FOPTs). However, it was not 
clear if the FOPTs observed in these models were just 
transients of DP. In context of non-equilibrium wetting 
process FOPTs are found in several ( 1 + 1 Vdimensional 
stochastic models with local interactions TO], contrary to 
equilibrium wetting process where phase transitions are 
not possible in 1— d systems having short-range inter- 
actions between interfaces and substrates. To the best of 
our knowledge FOPTs in deterministic, chaotic, extended 
systems with short range interactions are still lackin g, al - 
though it is known to exist for globally coupled maps |l2j | . 

In this Letter we introduce and study a single param- 
eter family of piece-wise linear chaotic maps which are 
diffusively coupled on a d-dimensional square lattice. A 
delay is introduced dynamically between sub-lattices by 
updating them alternately. One of our interests would 
be to find if, starting from a random initial condition, 



these sub-lattices synchronize at later times. The answer 
turns out to be 'no', for both very high and low diffu- 
sion strengths e. However for intermediate e synchro- 
nization occurs with the suppression of spatio-temporal 
chaos. This synchronized phase is an unique absorbing 
state of the system and for piece- wise linear maps (PLMs) 
the phase boundary is identical with the boundary for a 
stable fixed point. 

For e = 0, the system visits the whole phase-space 
in time. The phase-space volume decreases as e is in- 
creased and at a critical strength eA, the phase-space 
suddenly shrinks to a point which corresponds to the 
fixed point of the primitive map. This fixed point is sta- 
ble until e = £b > sa- Further increase in e results in 
de-synchronization where the phase space splits into two 
ergodic regions about the collective bifurcation points. It 
may be argued that for PLMs the width of the bifurcation 
is independent of e and thus synchronization error is dis- 
continuous at the critical point. However for non-linear 
differentiable maps the width of the bifurcation vanishes 
as ^/e — Cc which results in a continuous transition. We 
argue that this continuous transition is in a different uni- 
versality class than that of KPZ and DP with critical 
exponents /3 = 1/2, j3/v = and 7 = 1 in 1-dimension. 

The Model : Consider a d-dimensional hyper-cubic lat- 
tice C of coupled identical maps /(m, zj), where is a 
real variable at site i = (ii, 12 . . . id) with ik varying from 
1 to L. We define the even and odd sub-lattices {C^ and 
C° respectively) as C^'" — {i : J2k^k = even, odd}, and 
denote (y^) as the variable of {Cf)- Starting from a 
random initial configuration, {a;^} and {y^} are updated 
alternately as 
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where A/j is a set of 2d nearest neighbors of i and e is the 
coupling strength, can be seen as a diffusion constant. 
Equivalently, in the first half unit of time {2;^} are up- 
dated while {j/j} are at rest and the reverse happens in 
second half. We will see later that the delay which is in- 
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troduced dynamically between sub-lattices is responsible 
for a complete synchronization of the system. Note that 
periodic boundary configuration in all d dimensions are 
used throughout. 

Synchronization occurs when the difference between 
and its neighbors vanish at all sites as t ^ oo. Thus, the 
order parameter can be defined as (j) — (</>*) where 
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and steady state average is taken over time and realiza- 
tions. Obviously vanishes in the synchronized phase 
and in the unsynchronized phase > 0. A trivial syn- 
chronized phase would correspond to the stable fixed 
point of the CLM, i e, {^^ = z*}. It is easy to see from 
Eq. ^ that z* — f{m,z*). For chaotic CLMs without 
delay the largest Lyapunov exponent is independent of 
e and is the same as the Lyapunov exponent of primi- 
tive map which is positive. Hence, a fixed point solution 
{zj = z*} is unstable for any e. With a delay, however, it 
can become negative in a region es < e < e^i- For PLMs 
CA.B corresponds to the boundary of linearly stable re- 
gion which can be calculated as follows. 

Let us take the initial state to be close to the fixed- 
point, i.e., = z* + (5xj and = z* + (5yj. In Fourier- 
space the evolution of Sx-j' and (5yj reads as. 
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where Rr — 4-, 
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H = f{m,z*) 



and jl — — e). Let E± denote the eigenvalues of the 
matrix defined in Eq. ||2Jl. From the stability require- 
ments \Max{n{E+),n{E^))\ < 1, 

we find e^i = 1 — and = (/i + l)/(2/i) which are 
drawn in the inset of Fig. ^as a phase boundary for the 
synchronized phase. For simplicity, it is assumed here 
that primitive maps have only single nonzero fixed point 
z*. One can further generalize it to maps with more fixed 
points. 

To find out the behavior of cj) close to these transitions 
we first restrict ourselves to one dimension and study a 
specific single parameter family of maps : 
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_ j mz/(m — 1) z < a 
m(l — z) z > a 



(4) 



where o = 1 — 1/m. This piece- wise linear mapping of 
[0, 1] onto itself is everywhere expanding for m > 1, and 
thus chaotic, with an invariant density uniform on [0, 1]. 
A particular example of this family with m = 2 is known 
as tent map. Note that the fixed point is z* = m/{jn+l). 

Synchronization : Let us first discuss the transition 
from the unsynchronized phase A (inset of Fig. ^ to 




FIG. 1: Numerically obtained (p is plotted against e for tent 
map {fi = —2). The inset shows the phase diagram in e-fi 
plane. 
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FIG. 2: This figure shows how the phase space changes in zi- 
Z2 plane as e — > ea; equivalently 5^0. For large S, the phase 
space is identical for two different initial configurations and 
82. However, dynamically different shapes are generated as 
5^0. The symbol '*' represents the fixed point 2* — .6. 



the synchronized phase. Close to the transition we take 
e = eA — S and find that the system become multi-stable 
as (5 — > 0, i.e., there are large number steady states and 
one out of them is chosen by the system depending on 
the initial configuration. Multi-stability has been dis- 
cussed earlier in the context of CMLs [13 , and systems 
of delayed differential equations |Ti|. It may be argued 
that the multi-stability is extensive, i.e., the number of 
attractors grow exponentially with the system size. Thus 
any statistical average has to be taken over large number 
of independent realizations, which resticts us to simmu- 
late large systems. We carried out numerical simmula- 
tions for L = 1024 and m = 2 {tent map) and find that 
(j> vanishes discontinuously at ea = 0.5 (see Fig. QJ. To 
confirm that it is a true first order transition, not just 
a transient effect, we monitor the phase space of every 
neighboring pair of co-ordinates as (5 — )■ 0. For example 
in Fig. [21we demonstrate how the phase space changes in 
Z1-Z2 plane. Every other pair of neighboring co-ordinates 
show similar changes. In practice, no noticable change is 
observed in the phase space when S < 10~* and then 
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suddenly the fixed point zi = z* = Z2 appears at 5 = 0. 
In other words, when (5 « 0, we have \zi — > for 
every realization and thus has a jump at the critical 
point CA- 

De- synchronization : The synchronized state persists 
up to e = e_B where phase space splits into two discon- 
nected ergodic regions. In the new phase B shown in the 
inset of Fig. ^ the system fluctuates about its collective 
bi-bifurcation points Xi = x* and yi = y* . Using Eq. 
we get 



X* ^ (l-e)/(x*)+e/(y*) 
y* - {l-e)f{y*) + efix*), 



(5) 



which can be solved for the family of maps defined in 
Eq. as = a± and corresponding y* = aip, with 

m^(2(5+l)±(2(5m+l) , r i-i J • 

"± = m-(2S+i)+m{iS+i) and (5 = 6 - £5. Dependmg on 
the initial configuration, different parts of the sub-lattices 
are then attracted to a+ or a_ with kink-like interfaces 
(inset (a) of Fig. separating them. It will be shown 
later that the width of such a kink w, diverges as 1/ VS as 
(5 — > 0. Thus, stable kinks can not be generated when d is 
0{L~^) and we have (f) = q;_|. — a_ . Clearly the jump in 
at the critical point es is A = 2/(m^ + m). Note that in 
the other limit, i.e., when L — > oo before 6 0, we have 
a slightly different (j). In this case the large systems will 
generate an average density of kinks, say p. Immediately, 
we have (j> = A(l — ^f^), where A is the area bounded 
by an even and an odd kink. Nevertheless the transition 
is discontinuous. 

The profile of a kink can be calculated as follows. Let 
us assume that it starts at site k with Zk = and 
Zk+i = a+. In steady state Zk+2 has two solutions : a- 

and a = m42^+i)+™(4^+i) ' obvious that the first 

solution, Zk^2 = ct- 1 does not generate a kink. Thus the 
steady state profile of a kink can be obtained as 



Zk+i+l 
Zk+i+2 
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where 9 = 2 tan^^ {y^ S / es) ■ Explicit form of Zk+i is 
lengthy to present here, rather we compare numerically 
obtained Zk+i with Eq. © in the inset (b) of Fig. 
|3| When site i is far from k and (5 w one can 
write an approximate solution as Zk+2i = ct- + and 

■ 2 c 

Zk+2i+i = a- - Ai, where A^ = 2m(m+i) ■ ^hc width 
of the kink is thus w = 2n, such that A„ = — a_. 
Clearly, w diverges as 1/VS. In Fig|31we have plotted w 
obtained from numerical simulations versus S in log-scale, 
which confirms this power-law. 

Higher Dimensions : The phase diagram for coupled 
PLMs in higher dimension is the same as shown in the 
inset of FigH Close to the transition point e^, phase A is 
mutistable and shows sudden shrink in phase space vol- 
ume. The desynchronization transition occurs at es . But 
unlike one dimension, the kink-type solutions in phase B 
are no more stable and hence jump in the orderparame- 
ter is A = 2/(m^ -I- to). Numerical simulations in 2 and 
3 dimensions confirms this finding. 
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FIG. 3: A typical steady state profile of odd (solid-line) 
and even (dashed line) sub-lattices are shown in inset (a). 
Inset (b) compares the kink profile obtained from simulations 
(points) with Eq. Q. In the main figure, width of the kink 
w is plotted against S in log scale for m — 2 and 2.5. The 
slope of the solid line is set to —1/2. 



Other Maps : How do non-linear maps behave when a 
delay is introduced between sub-lattices ? Clearly the 
phase boundaries for synchronized phase strongly de- 
pends on corresponding primitive maps. First let us con- 
sider the logistic map : f{x) = 4x{l — x). The fixed 
point of this map Zi — z* — 3/4 is linearly stable for 
1/2 < e < 3/4. However numerical simulation shows 
that the synchronization occurs at eA > 1/2. Close to 
the transition point phase A is multistable and hence we 
have a first order transition similar to that of the tent 
map. The de-synchronization is found to occur exactly 
at 7 where the linear-stability of the fixed point breaks 
down. Unlike the tent maps here we do not have fluc- 
tuations about the bifurcation points. But the kinks are 
unavoidable. Using Eqs. |(SJ) and © one can obtain 
the bifurcation points as Xi — a± and y^ — azp, where 

a± — 4{i+iS) ■ Clearly A — a+ — a_ vanishes as 

^/S indicating that the transition is continuous. Critical 
exponents of this continuous transition can be deduced 
using following arguments. 

First, note that fluctuation of order-parameter comes 
about from the variation in density of kinks which fluc- 
tuates about a mean density, say p. Density distribu- 
tion about p can be assumed to be normal, which gives 
(f) ~ VS{1 — ^^), where A is the area bounded by an 
odd and an even kink, and w is the width of the kink. 
Since <j) is independent of L and proportional to ^/S, we 
have critical exponents /3 — 1/2 and /3/i^ = 0. One can 
further define "susceptibility" as x = (('^^) — (<^)^)/-^- 
Now X ^ Sr'^/L, and hence 7=1. It is evident that the 
critical dimension dc = 1- 

We have also done numerical simulations on sine 
{f{x) = sin(7ra;)) and cubic {f{x) = V27x{l - x'^)/2) 
maps in one dimension, which show that the syn- 
chronization transition is discontinuous, whereas de- 
synchronization occurs continuously with the same criti- 
cal exponents as that of the logistic map. 

To get an insight why de-synchronization is continuous 
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FIG. 4: Phase diagram of de-synchronization transition is 
shown in e-h plane. The nature of transition changes from 
being discontinuous (dashed line) to continuous (solid line) 
at 6 = 1. In the inset, <j) and x versus 5 is plotted in log-scale 
and the solid lines are drawn with slope (3 = 1/2 and 7 = 1 
respectively for comparision. 

for certain class of maps we investigate coupled power- 
law maps: f{z) — 1 — \2z — Ij'', where 6=1 and 2 cor- 
responds to the tent and the logistic maps respectively. 
Thus, by tuning b it is possible to study how the na- 
ture of transition changes from being first to second or- 
der. Our numerical simulations suggests that for 6 > 1 
the de-synchronization transition is continuous and be- 
longs to the same universality class as that of the logistic 
map, whereas for 1/2 < 6 < 1 it is discontinuous. Note 
that there are two non zero fixed points for & < 1/2 and 
present analysis need to be modified for their discussion . 
Basically, synchronization error is the average difference 
between bifurcation points x* and y*, which vanishes at 
es where x* = 1/2 = y*. Just above e_B this difference is 
proportional to the difference between the left and right 
slope of f{z) at z = 1/2. For b < 1, since these maps 
are not differentiable at z = 1/2 the de-syncronization 
occurs discontinuously, whereas for b > 1 the maps are 
differentiable and thus we have a continuous transition. 
A few comments are in order. First, the synchronized 



phase discussed here is different from earlier studies, 
where CMLs remain chaotic in the synchronized phase 
for both, when they are coupled directly 0| or through 
noise . In our model the spatio-temporal chaos are sup- 
pressed and thus the synchronized phase is a completely 
ordered phase. The spatial correlations even persist in 
de-synchronized phase B. Second, one can define a differ- 
ent order parameter called "magnetization" M = (m*), 
where m' = sign{zl — z*). Then the desynchronyza- 
tion transition can be thought as a transition from fully 
ferromagnetic to fully antiferrmagnetic phase, which can 
be either first or second order. Note that the continu- 
ous transition obtained here is in a different universality 
class than that of zero temperature Ising transition in 
one dimension. 

Conclusion : In conclusion, we study a single parame- 
ter family of coupled chaotic maps by introducing a de- 
lay between sub-lattices. We show that these systems 
undergo a first order transition to a synchronized phase 
when the coupling parameter is varied. The transition 
occurs as the system enters from a multiply stable re- 
gion to a single 'collective fixed point' in phase space. A 
second transition occurs as this fixed point become un- 
stable and the phase space breaks up into two ergodic 
regions about the collective bifurcation points. Our an- 
alytical results show that the de-synchronization transi- 
tion is discontinuous in one and higher dimensions when 
the primitive maps are not differentiable. Differentiable 
maps also show a first-order synchronization transition, 
whereas the de-synchronization occurs continuously with 
critical exponents (3 — 1/2, 7=1, and = 0. We ar- 
gue that the de-synchronization transition discussed here 
can also be thought as a transition from a fully ferromag- 
netic to a fully anti-ferromagnetic phase. These results 
establish a true non-equilibrium first order transition in 
one dimensional dynamical systems with short range in- 
teractions. 
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